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Abstract 

For a system of fermions with a three-body contact interaction the second-order contri- 
butions to the energy per particle E{kf) are calculated exactly. The three-particle scatter- 
ing amplitude in the medium is derived in closed analytical form from the corresponding 
two-loop rescattering diagram. We compare the (genuine) second-order three-body contri- 
CS| ' bution to E{kf) ~ k^P with the second-order term due to the density-dependent effective 



/ 
two-body interaction, and find that the latter term dominates. The results of the present 

study are of interest for nuclear many-body calculations where chiral three-nucleon forces 

are treated beyond leading order via a density-dependent effective two-body interaction. 



PACS: 12.38.Bx, 21.30.Fe, 24.10.Cn 



1 Introduction and summary 

1^ ■ Recent advances in the formulation and construction of (low-momentum) nuclear interactions 
in chiral effective field theory have unambiguously revealed the important role played by three- 
nucleon forces [I1I21IS]- Three-body forces turn out to be an indispensable ingredient in accurate 
calculations of few-nucleon systems ^ as well as for the structure of light nuclei [S]. In chiral 
cn ' perturbation theory the three-nucleon interaction can be constructed systematically and con- 
^ ■ sistently together with the nucleon-nucleon potential [HE]. At leading order it consists of a 
\0 , zero-range contact-term, a mid-range Ivr-exchange component and a long-range 27r-exchange 
cn component, where the parameters of the latter component occur also in the (sub leading) 2tt- 
exchange NN-potential. The calculation of the subleading chiral three-nucleon force, built up 
from many pion-loop diagrams etc., has been completed recently in ref. [B] and applications to 
few-nucleon systems are underway. 
k> . Furthermore, it has been demonstrated that by employing low-momentum two-body interac- 

j_j ' tions (instead of traditional hard-core NN-potentials) the nuclear many-body problem becomes 
- . .' significantly more perturbative |7]. This desired simplification is accompanied by a prominent 
role of the three-nucleon interaction, such that its inclusion is essential in order to achieve 
(reasonable) saturation of nuclear matter already at the Hartree-Fock levelo The combined 
repulsive three-body effects counterbalance with increasing density the purely attractive con- 
tributions provided by the low-momentum two-body interactions alone. Improved calculations 
of nuclear matter which aim at reproducing the empirical saturation point, £'0 — — 16MeV, 
Po — 0.16 fm~^, still have to treat second-order (and even higher-order) corrections which arise 
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^At this point it should noted that three-nucleon forces are unavoidable, independently of whether one uses 
low-momentum or conventional NN-potentials for nuclear matter calculations. Hard-core NN-potentials are 
simply not usable in perturbative many-body calculations and one has to treat them at least at the Brueckner- 
Hartree-Fock level. 
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Figure 1: Left: Three-body contact coupling modelled by heavy scalar-isoscalar boson exchange. 
Right: Closed three- loop diagrams representing the energy density linear in C3. 



in many-body perturbation theory from the low-momentum two- and three-nucleon interac- 
tions [S|. An approximate treatment at second order (and beyond) is commonly pursued by 
mapping the three-nucleon force onto a density-dependent effective two-body interaction. The 
detailed form of the density- dependent effective NN-interaction as it results from the leading 
order chiral three-nucleon interaction has been worked in ref. p], in particular with regard to 
its implementation into nuclear structure calculations. 

The present paper aims to overcome this common approximation by performing an exact 
second-order calculation for the simplest three-body interaction, namely for the zero-range 
contact interaction. Our paper is organized as follows. In section 2 we recapitulate the first 
order calculation of the energy per particle E{kf) ~ k^ and outline different methods to compute 
the spin-isospin weight factors of closed three-body diagrams. In section 3 the two-loop diagram 
describing the three-particle rescattering in the medium is evaluated in detail. Analytical 
expressions are derived for the corresponding vacuum term Bq, and for the medium corrections 
Bi and B2 incorporating Pauli-blo eking effects due to one and two particles, respectively. The 
ultraviolet divergence in the vacuum loop Bq requires the introduction of a (Galilei-invariant) 
three-body counterterm proportional to fourth powers of momenta. The analytical results 
for Bqi2 are then used in section 4 to compute the second-order (i.e. five-loop) contribution 
to the energy per particle E{kf) ~ fc}°. By comparing this novel result with the second- 
order term provided by the density-dependent effective two-body interaction, one deduces from 
the numerical prefactors that the latter term actually dominates (by about a factor 2). The 
appendix contains two remarkable reduction formulas for special integrals over the product of 
three Fermi spheres. 

The present result obtained in an exact second-order calculation with a simple three-body 
contact interaction may serve to support the approximate treatment of chiral three-nucleon 
forces in nuclear matter calculations [71 18]. Nevertheless, complete second-order calculations 
with the finite-range chiral three-nucleon forces should be attempted in the future. 



2 First order calculation and weight factors 

We start out with reproducing the first order calculation of the energy per particle E{kf) from a 
three-body contact interaction. This allows us to fix the notation for the coupling constant and 
to discuss alternative methods for computing the spin-isospin weight factors of closed three-body 
diagrams. The involved combinatorics becomes most transparent if one models the contact- 
vertex C3 by the exchange of two heavy scalar-isoscalar bosons. The topologically distinct 





Figure 2: First and second order diagrams generated by a three-body contact interaction. 

diagrams obtained by closing all three nucleon lines are shown in Fig. 1 (together with symmetry 
factors 1/2). The resulting spin-isospin weight factor in this approach is 4^/2—4^—4^/2+4 = 12, 
and hence the first order contribution to the energy per particle reads: 
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with the density given by p 
sphere. For orientation and comparison we note that the coupling constant of the contact- 
term in the chiral three-nucleon force [T] is parameterized as C3 = ce/ f^^x- '^^^ reduced 
coupling constant Fs = —C^/12ti'^ is particularly useful to write down compactly the second 
order contributions (see eqs. (14,16)). 

An alternative derivation of the spin-isospin weight factors starts from a contact-vertex that 
additionally includes the three-particle antisymmetrization operator: 

A=Y. sign(a)P,[o-T] = (1 - Pi2)(l - P13 - ^23) • (2) 

Here, Vij = (1 + ctj ■ (Tj)(1 + Tj ■ tj)/A is the exchange operator for particles i and j, operating 
in their respective spin and isospin spaces. For the first order three-body diagram (shown in 
the left part of Fig. 2) a trace has to be taken over the spin- and isospin degrees of freedom of 
all three nucleons, leading to the result: 

^ Tr^ = ^ 64 f 1 - 1 - i - i + 1 + lU ^ 64 ■ ^ = I2C3 . (3) 

6 2V44416 16/2 8 ^ ^ ' 

Here, the symmetry factor 1/6 accounts for the permutations of the three (indistinguishable) 
closed nucleon lines and we have displayed separately the contributions from all six permuta- 
tions a G 5*3. The effective coupling constant is now 3C3, since in this completely symmetric 
formulation of the contact- vertex each of the three incoming particles must be joined to the cen- 
tral vertex of the boson-exchange picture in Fig. 1. More interesting is the spin-isospin weight 
factor of the second order diagram shown in the right part of Fig. 2. Taking into account its 
symmetry factor 1/6^, one gets: 

i^^Tr^2 ^ Cl^^ g^ ^ 3g^2 ^ 12C3 • 3C3 , (4) 

6^ 4 

where the relation J\? = 6 A has been used. The factorization at the end of eq.(4) shows that 
the three-particle rescattering process carries the relative weight factor 3C3. 

We finish this section by considering the density-dependent effective two-body interaction 
which results from closing one nucleon line of the three-body contact- vertex (see Fig. 3). The 




Figure 3: Density- dependent effective 2-body interaction obtained by closing one nucleon line. 
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Figure 4: Two-loop diagrams describing the three-particle rescattering in the medium. 

method of employing the boson exchange picture gives for the two-body coupling strength: 



5C,{p) = C, 
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where the factor 6 = 4+4-2 — 4 — 2 emerges from the set of 9 topologically distinct diagrams. Of 
course, the same result is obtained by including the (two-particle) antisymmetrization operator: 
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where tra sums over the spin and isospin degrees of freedom of the third nucleon. Having the 
expression for the effective two-body coupling, 5Co{p) = Cs/cj/vr^, one can immediately give 
the result for its second order contribution to the energy per particle (see eq.(16)). 



3 Three-particle rescattering in the medium 

The calculation of the second-order (five-loop) diagram in Fig. 2 proceeds via the three-particle 
scattering amplitude. The two-loop diagrams describing the three-particle rescattering in the 
medium are shown in Fig. 4. By including the topologically distinct diagrams with crossings 
of internal lines the crucial factor 3 deduced in eq.(4) gets accounted for. Each internal line 
introduces a particle propagator, ^[1 — 6{kf — \lj\)]/{loj — l^/2M + ie), and after performing 
the energy integrals via residue calculus the expression for the in-medium scattering amplitude 
takes the form: 



B 



d%d% 
(27r)6 



3C,,M 



I? + 3/1/4 - H/Q - It 



-e{kf-\p + h-h/2\) 
eikj~\p-h-k/2\) 



(7) 



The chosen assignment of intermediate momenta: p + li — hf^, p + h, P — h — ^2/2, with 
P = {Pi +P2 + P3)/3, has the advantage that interference terms of the loop momenta li^2 among 
themselves and with the external momenta pj are absent in the energy denominator. The 
external momenta are to be taken from the region inside the Fermi sphere \pj\ < kf, hence 
the dimensionless variable s = \p\/kf satisfies the condition < s < 1. The other kinematical 
quantity H appearing in the energy denominator is the Galileian invariant: 

H = {pi- P2f + {pi - nf + (P2 - nf < dkj . (8) 

The maximum value 9/c? of H is reached in the configuration where ^1,2,3 point from the center to 
the vertices of an equilateral triangle of side- length \/3kf. The in-medium scattering amplitude 
B is manifestly real- valued, since Pauli-blocking and energy conservation forbid any imaginary 
part: 3k} < \p + k\^ + \p + h - h/^l^ + \p- h - k/'^? = (2? + 3/|/2 - H/?>) +p^+pi + p^ = 
Pi + pf + Pf < 3/:;?, where the term in brackets vanishes on-shell. 

In the next step one expands the product of the three (1 — ^)-factors in eq.(7). This 
rearrangement gives a vacuum part Bq with no ^-factor, the sum Bi of three equal terms 
with one ^-factor, and the sum B2 of three equal terms with two ^-factors. The equality of 
the three summands is obvious from the structure of the two-loop diagram and can be shown 
explicitly by making the substitution li = A/3/3/2. The intermediate momenta are then linear 
combinations of I2 and ^3 with coefficients given by rotations about angles ±27r/3, and the 
energy denominator 3{l2 + /3^)/4 — H/6 is clearly invariant under these discrete rotations. At 
last, there is a contribution B^ involving three ^-factors. Fortunately, one does not need to 
evaluate the complicated term Re-B3, because its contribution to the energy per particle E{kf) 
vanishes identically. The corresponding integral over six Fermi spheres is symmetric under the 
interchange of the external and internal momenta, except for the energy denominator which 
changes its sign. To be precise, the denominator in eq.(7) is interpreted in this argument as a 
principal value. 

Next, we have to evaluate the two-loop integrals for i?o, -Bi and B2. The divergent vacuum 
loop Bq is treated by dimensional regularization. The three-dimensional integrals in eq.(7) are 
continued to d dimensions by the rule: {2ti)~^ ^ d^lj — )■ X^~'^{27i)~'^Jd'^lj, which introduces a 
scale A in order to preserve mass dimension of the loop integrals. The divergent behavior of 
the vacuum loop Bq shows up through an Euler Gamma-function r(l — c?) and after expanding 
around d = 3 one gets: 
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A three-body counterterm proportional to H^ (i.e. fourth power of momenta) is needed to 
renormalize the three-particle scattering amplitude in vacuum. This is different to the case of 
two-body scattering, where the vacuum divergence (not visible in dimensional regularization) 
can be absorbed on the scattering length [TUIIII]. It is convenient to introduce via the relation 
H = dk'jh a second dimensionless variable h, which satisfies together with s = \p\/kf the 
constraint s"^ + h < 1. Using this variable the real part of the renormalized vacuum loop reads: 

.,p„^ 3V3CsMki „rl, , kf , J 

ReBt^ = "^^/^^ ^ h^ (3 - In 3) - 2 In -^ - In /i + ct A , 10 

[Air)'^ 12 A J 

with ct(A) a parameter for the scale-dependent counterterm. 

In order to evaluate the two-loop integral Bi we select the step-function 9{kf — |p + /2|)- 
In dimensional regularization the /i-integral is finite and proportional to v9l2 — 2H. The 



remaining integral over a shifted Fermi sphere of radius kf leads to the following result for the 
real part of Bi. 
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For the evaluation of i?2 we choose the product 9{kf — \p+li—l2/2\) 6{kf — \'p—li—l2l2\) of two 
step-functions. The integral over Zi-space leads to the familiar hole-hole bubble [IHIIII] involving 
several logarithms. The condition for it not to vanish is |p — /2/2I < kf and consequently the 
integration region in /2-space is a shifted Fermi sphere of radius 2k f. In order to get a more 
concise representation of B2 we introduce the following auxiliary function: 

1— N 

A(Q, A^) = 2JQ arctan -^= , for Q > , 
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The final result for the real-part of B2 written in terms of the variables s and h has the form: 
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Note that the symmetrization prescription -|-(s — ?■ — s) applies to all terms in the curly brackets, 
also the first few which are already even functions of s. We remark again that individual 
components Bj of the three-particle in-medium scattering amplitude possess an imaginary 
part, but these add up to zero: Im {Bq + Bi + B2 + B3) = 0. 

4 Results: Energy per particle at five-loop order 

Having available the analytical expressions for Rei?o,i,2 in eqs. (10-13) it is straightforward to 
compute the second order contribution to the energy per particle E{kf). The pertinent integral 
over three Fermi spheres (suitably parameterized by three radii and three angles) can be solved 
in closed form for the polynomial H^. Extracting the coefficient of this piece, the expression 
for the energy per particle at five-loop order (see right diagram in Fig. 2) reads: 

^(M = ^r|^^)°|v^ln^ + Co + Ci + C2}<0, (14) 
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Figure 5: Left: Density-dependent effective two-body interaction to second order. Right: 
Anomalous three-body diagram which vanishes at zero temperature. 

where the counterterm ct(A) has been absorbed into the logarithm ln(fcj/Ao). The values of 
the numerical constants Co, Ci and (2 as obtained by integrating over three unit spheres are: 

Co = -1.425, Ci = -5.653, C2 = -4.354 ± 0.014 , (15) 

where the errors of Co and Ci he beyond the digits given. The result in eq.(14) is to be compared 
with the second-order contribution from the density- dependent effective two-body interaction 
as represented by the left diagram in Fig. 5. Converting 6Co{p) into a scattering length a{p) = 
C^Mk^/AiT^ and respecting the spin-isospin degeneracy factor 4 — 1 = 3, one finds from the 
familiar low-density expansion [TU| [TT|: 

E{kf) = —TlMkf(ll-2\n2)>0. (16) 

The numerical factors multiplying T^Mkl^ in eqs. (16,14) have the values 54(11 — 21n2)/35 = 
14.83 and 377r(Co + Ci + C2)/175 = —7.59. Their comparison reveals that the repulsive second 
order contribution from the density- dependent effective two-body interaction is dominant, but 
gets reduced to about half of its size by the (genuine) second order three-body contribution. 
The additional logarithmic term — 1.151n(Ao//c/) does not change this balance in a significant 
way, assuming that the scale Aq is of natural size. For the sake of completeness we note that 
at five-loop order one can additionally construct the anomalous three-body diagram shown in 
the right part of Fig. 5. It involves in the central loop the squared in-medium propagator and 
it vanishes (at zero temperature) by reason of the identity: 9{kf — \l |)[1 — 9{kf — \l |)] = 0. 

The expressions in eqs. (14, 16) for the energy per particle E{kf) ~ r|M/c}'' at five-loop order 
constitute the main results of the present work. These analytical results demonstrate that by 
treating a three-body (contact) interaction through the density- dependent two-body interaction 
only [3 |8] second order effects in many-body perturbation get considerably overestimated. 
Clearly, after having made this observation in a special case one should attempt complete 
second-order calculations with the finite-range chiral three-nucleon forces in the future. 

Appendix: Reduction formulas for integrals over three Fermi spheres 

In this appendix we present two remarkable reduction formulas for integrals over the product 
of three Fermi spheres. In the first case we assume a dependence of the integrand on the 
symmetric variable s = \pi + p2 + 'P3\/3kf. The following reduction formula holds: 

d^Pid^P2d^P3 ^f \Pi+P2+P3\ \ ^ 9fc| n 
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Figure 6: Weight functions w(s)/8 and xix). The area under both curves is equal to 280/243. 

This peculiar result for wi^s) has been obtained with the help of Fourier transformation tech- 
niques and eq.(17) has been checked numerically for many examples of Fi^s). 

In the second case we assume a dependence of the integrand on the combination s^ + /i = 
[pi +p| +p|)/3/i;? and derive the following reduction formula: 
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with the (non-smooth) weighting function 
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This intricate expression for x{^) has been obtained by inverting the order of integrations over 
the three involved radial coordinates p'j/Sk'j, while fixing their sum to x. Numerical checks of 
eq.(19) have been performed for many examples of G{x) as well. 

Fig. 6 shows the two weight functions w{s)/8 and x{^) next to each other. The area under 
both (bell-shaped) curves is equal to 280/243 = 1.152. A quick analysis gives that w{s)/8 
reaches its maximum value of 2.65 at s = 0.393, while x(a;) reaches its maximum value of 2.92 
at a; = 0.621. The behavior at the kinematical endpoints is opposite: w{s) = 350s^, s — ?■ and 
w{s) = 1890(1 - s^, s -)■ 1 versus xi^) = Gt^V^x'^^'^, a; -)■ and x{x) = 105(1 - a;)^/2, a; -)> 1. 
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